Discovery of Molecular Hydrogen in White Dwarf Atmospheres 
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ABSTRACT 

With the Cosmic Origins Spectrograph onboard the Hubble Space Telescope, 
we have detected molecular hydrogen in the atmospheres of three white dwarfs 
with effective temperatures below 14,000 K, G29-38, CD 133 and CD 31. This 
discovery provides new independent constraints on the stellar temperature and 
surface gravity of white dwarfs. 

Subject headings: white dwarfs, atmospheres 



1. INTRODUCTION 

An important discovery from the International Ultraviolet Explorer (lUE) was two broad 
absorption features at 1400 A and 1600 A i n hydrogen-dominated (DA) white dwa rfs cooler 



than 20,000 K and 13,500 K, respectivel v flGreenstein fc Okel ll979 



Wegnei]ll982h. Subse- 



quently, these features were explained by iKoester et al.l (1l985[ l and iNelan fc Wegnerl (1l985[ l 
as Ly a satellite lines from collisions among H-H (1600 A) and H-H+ (1400 A). These quasi- 
molecular abs orption lines help to co nstrain stellar temperature and surface gravity for DA 
white dwarfs JAUard fc Koestei]ll992l ). 029-38, GD 133 and GD 31, the targets described in 
this letter, a ll show these quasi-mole c ular a bsorption features when observed with the WE 



Holm et al.l ( 119851 ): iKepler &: NelanI (jl993[ ). lUE archive]. Here, we report first detections 



of true molecular hydrogen in these three white dwarfs, by employing the Cosmic Origins 
Spectrograph (COS) on the Hubble Space Telescope (HST). These three stellar atmospheres 
are among the hottest stellar environments where photospheric molecular hydrogen has ever 
been detected. 
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DATA 



G29-38 and GD 133 were observed as part of program 12290 during HST cycle 1 
which focuses on using externally pollut ed white dwarfs to assess the v olatile abundances 



of acc r eted extrasolar plari e tesimals [see lJural (|2003[): IZuckerman et al.l (120071 ): iKlein et al. 
J2010h ; bufour et all (EolJ ; ljura et all J2OI2I ): ICansicke et al.l J2012h and references within] 



G29-38 and GD 133 were chosen because they show a high degree of atmospheric poUu- 
tion from optical studies (IKoester et al.l 19971 l2005l) and displa.y exc ess infrared radiation 



(IZuckerman fc BecklinI 119871 : iReach et al. 



gram to study Hyades white dwarfs (IZuckerman et al. 



12005 



2009 



Jura et al.ll2007l ). As part of our pro- 



2OI3I ). ultraviolet data for GD 31 were 



retrieved from the i/ST archive of the SNAPSHOT program 12169 (PI: B. Gansicke). 

For G29-38 and GD 133, the C OS set-up was s imilar to that employed for GD 40 
and G241-6, as previously reported in I Jura et al.l ( 120121 ). The G130M grating was used with 
central wavelength 1300 A and wavelength coverage of 1142 -1288 A (strip B) and 1298 -1443 
A (strip A). The spectral resolution was ~ 20,000. Total exposure times, from combining 4 
separate exposures for G29-38 and 5 for GD 133, were 9032 sec and 13,460 sec, respectively. 
The signal-to-noise ratio is at least 8 for strip B and 15 for strip A for both stars. The raw 
data were processed under the pipeline CALCOS 2.18 .5. Because G29-38 is a ZZ Ceti star 
with multiple pulsation modes ( iThompson et al.ll2008l ). which lead to considerable changes 
in flux levels and continuum shape in each individual exposure, CALCOS processing cannot 
fully correct for these factors. Instead, we fitted the continuum of each CALCOS-extracted 
spectrum with a low-order spline3 polynomial and combined the normalized spectra with 
IRAF. We did not find any noticeable differences in line strengt h in each exposure. GD 
133 is also a ZZ Ceti but the pulsation amplitude is much lower (ISilvotti et al.ll2006l ) and 
we simply adopt the pipeline reduced spectra. For GD 31, th e G13 0M configuration has 
a central wavelength of 1291 A, as described in iGansicke et al.l (120121 ). The total exposure 
time was 1200 sec and the signal-to-noise ratio is ~ 7 for strip A. We also used the output 
from the CALCOS pipeline for GD 31. 

Apart from th e absorption lin es from heavy elements in G29-38 and GD 133, which 
will be discussed in IXu et al.l (l2013l ). we see numerous quasi-periodic features, as illustrated 
in Figure 1. They are ubiquitous between 1310 A and 1443 A, almost the entire strip A. 
Clearly, these features are real and not instrument al artifacts, such as fixed pattern noise, for 
the following reasons. Unlike GD 40 and G241-6 (jJura et al.ll2012l ). two other white dwarfs 
observed in the same program with the same instrumental set-up, these features were only 
present in the spectra of G29-38 and GD 133. Furthermore, they are at the same wavelength 
in each separate expo sure. Although initi ally puzzled, after comparing with the ultraviolet 
spectrum of the Sun (ISandlin et al.l 1 19861 ). we realized the absorption wavelengths exactly 



- 3 - 



match those of molecular hydrogen. All the recognizable features correspond to Lyman 
band tran sitions from v"= 2, 3 , 4, 5 to z/' = and Werner band transitions from v"=2, 3 
to z/' = ( lAbgrall et al.lll993al Jbl). It turns out that we have accidentally found H2 in white 
dwarf atmospheres! We present the strongest Lyman band H2 lines in Table [T| they are all 
appreciably stronger than the Werner bands in our observed wavelength interval. 

We serendipitously identified molecular hydrogen in a third star, GD 31 , which may be 
a high-mass escaping member of the Hyades cluster (jZuckerman et al.l 120131 ). The data are 
noisier but clearly four absorption features are seen; these correspond to the strongest H2 
Lyman band lines and their blends in G29-38, at the correct wavelengths, as presented in 
Table [Hand Figure 1. 



DISCUSSION 



According to the Saha equation, molecular hydrogen is most concentrated in low tem- 
perature, high density environments. Previously, ro- vibrational lines from H2 have been 



detected in cool stellar atmospheres in the infrared (ISpinrad 



19641 ) and in the Sun's ultravio- 



let emission spectrum ([Jordan et al.lll978l : ISandlin et al.lll986| ). Our ultraviolet detection of 
photospheric H2 introduces a ne w aspect to stel lar physics. With stellar parameters of T* = 
11,820 K, log g=8.40 for G29-38 JXu et al.ll2013[ ) and T,=12,121 K, l og g = 8 .005 for GD 133 



(IKoester et al.[[2009[ ). we computed white dwarf model atmospheres (lKoester|[2010[ ) and cal- 
culated the number density of molecular hydrogen relative to atomic hydrogen, n(H2)/n(H), 
as shown in Figure 2. At maximum, molecular hydrogen is still 10~^'^ and 10~^'^ less than 
the amount of atomic hydrogen in G29-38 and GD 133, respectively. These computed ratios 
are comparable to, but at the higher en d of the concentration of trace elements in heavily 
polluted white dwarfs (I Jura et al.|[2012[ ). Because molecular hydrogen is distributed over a 
large number of rq-vibr ational levels, each individual line is relatively weak. As discussed in 
Zuckerman et al.l (120131 ) . H2 is used to resolve a major temperature puzzle for GD 31. Due to 
its high gravity, the equivalent widths (EWs) of 4 detected H2 lines in GD 31 (T=13,700 K, 
log g=8.67) are comparable to those in GD 133; there is a substantial amount of molecular 
hydrogen in the atmosphere of GD 31. 

As s hown in Figure 1, we computed the model spectra for G29-38 and G D 133 following 
Koester (2010) with all the line dat a obta i ned fr om the Kurucz webpag^ (Kurucz 1995 ) 
and the H2 partition function from llrwini (jl98ll ). In a statistical sense, the model well 
represents the data and reproduces most molecular hydrogen features. However, the fit to 
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individual lines is rather poor due to the lack of accurate broadening parameters. The Z 
Astrophysical Plasma Properties (ZAPP) Collaboration is actively involved in c reating white 



dwar f photospheres in the lab and determining line broadening parameters (iFalcon et al. 



20121 ). Our understanding of the COS ultraviolet observations can be substantially improved 



with accurate physical parameters for the H2 lines. 

Another aspect of this discovery is that by constraining the abundance of the HD 
molecule, we can place a limit on the D/H ratio. It is generally believed that all deu- 
terium is prirnordia l from the Big Bang nucleosynthesis and destroyed in stellar interiors 



(lEpstein et al.lll976[ ). Any amount, if present, in the white dwarf atmosphere must come 
from some external source, likely relic planetesimals. The values are not very constraining for 
the targets presented here. But if cooler stars can be observed with COS with a sufficiently 
high signal-to-noise ratio, then more meaningful upper limits or actual detections of the HD 
molecule in white dwarf atmospheres may be anticipated. 

In white dwarfs cooler than 12,000 K, H2 is present but usually there is not enough 
ultraviolet flux to make observations in a time-efficient manner. The environment in hotter 
stars is typically more hostile for molecular hydrogen but high pressure may still enable 
detection of H2. 



4. CONCLUSIONS 

With the HST, molecular hydrogen was detected for the flrst time in white dwarf at- 
mospheres. The three stars, G29-38, GD 133 and GD 31, have temperatures between 11,800 
K and 13,700 K. H2 can be used as an independent constraint to white dwarf atmospheric 
conditions. This opens a door to many future explorations. 

Support for program # 12290 was provided by NASA through a grant from the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. We also utilized observations 
obtained [from the Data Archive] at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under NASA contract 
NAS 5-26555. These observations are associated with program # 12169. This work also has 
been partly supported by NSF grants to UCLA to study polluted white dwarfs. 
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Fig 1. Portion of the Lyman band absorptions of molecular hydrogen in G29-38, GD 
133 and GD 31 at the red wing of the H-H+ quasi-molecular feature. The spectra are 
smoothed with a 5 point boxcar and shifted by 39 km s~^, 58 km s~^ and 89 k m s~^ fo r 



G29-38, GD 133 and GD 31 respectively, to be in the heliocentric reference frame (IXu et al. 



20131 : IZuckerman et al.ll2013l ). For clarity, the spectrum for GD 133 is offset by 0.5 x 10 



-14 



erg s~^ cm~^ and for GD 31 by 1.5 x 10~^^ erg s"^ cm^^ A~^. The red dashed lines 
are our computed model spectra. Statistically, the models reproduce the data but the fit 
to individual lines is not ideal. For G29-38 and GD 133, the entire spectrum in strip A 
(1298-1443 A) shows numerous ro-vibrational lines of H2. The strongest absorption features 
around 1345 A, 1357 A, 1366 A and 1369 A are also seen in GD 31. 
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Fig 2. Ratio of number densities of molecular hydrogen compared to atomic hydrogen 
at different optical depths in G29-38 and GD 133. For G29-38, this ratio peaks at log 
n(H2)/n(H) = -4.7, where r = 0.04, T = 9900 K and p = 3.1 x 10"^ g cm'^. For GD 133, 
the maximum of log n(H2)/n(H) is -5.1, which occurs at r = 0.05, T=10,200 K and p = 1.5 
X 10"'' g cm~^. 
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Table 1: The Strongest Identified Lyman Band H2 Lines 



A 


A 


Transition'' 


G29-38 


GD 133 


GD 31 




(A) 




EW (mA) 


EW (mA) 


EW (mA) 


74344.49 


1345.1 


0-4 R(5) 


170" 


107" 


87" 


74339.87 


1345.2 


0-3 P(ll) 


170" 


107" 


87" 


73936.94 


1352.5 


0-4 P(5) 


76 


66 




73719.80 


1356.5 


0-4 R(7) 


121" 


103" 


63" 


73714.96 


1356.6 


0-3 P(12) 


121" 


103" 


63" 


73185.32 


1366.4 


0-4 P(7) 


180 


70 


59 


73064.05 


1368.7 


0-3 P(13) 


45 


35 


35 


71225.99 


1404.0 


0-4 P(ll) 


76 


45 




71187.04 


1404.8 


0-5 R(5) 


97 


83 




70992.76 


1408.6 


0-3 P(16) 


112" 


90" 




70889.48 


1410.6 


0-4 R(13) 


112" 


90" 




70076.46 


1427.0 


0-4 P(13) 


75 


53 





" From lAbgrall et al.l fll993af ). 

We follow the notational conventions in iMeyer et al.l (2001). 
" This line is blended and the reported number is the total EW for both lines together. 
Note. For G29-38 and GD 133, the EW uncertainty is dominated by the choice of continuum 
interval because the whole region shows numerous H2 lines and there is no clean continuum 
as shown in Figure 1. Typically, the EW uncertainty is within 10%. For GD 31, the data 
are much noisier and the EW uncertainty is typically 25 %. 



